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(57) Abstract: An MRI imaging system indudes at least one processor and a plurality of coils to acquire a plurality of k-space 
2 samples of a taiget to image. The system includes a machine-rBadable media comprising instrucdons which, when executed by die 

processor, result in determining a plurality of different regularization matrices for a plurality of different regions of an image of the 
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A METHOD AND APPARATUS TO IMPROVE AN MRI IMAGE 

RELATED APPLICATION 

5 This application claiins the benefit of U.S. Provisional Patent Application 

- No. 60/348,005 filed October 19, 2001. 

TECHNICAL FIELD 
The disclosure relates generally to Magnetic Resonance Imaging (MRI), and 
10 more particularly to arti&ct cancellation in MR generated images. 

BACKGROUND 

Magnetic Resonance Imagiiig (MRI) is an imaging technique based in part 
on the absotpticm and emission of energy in the radio fiequency range. To obtain 
the necessary magnetic resonance (MR) images, apatient (or other target) is placed 

IS in a magnetic resonance scanner. The scanner provides a magnetic field that causes 
magnetic moments in the patient or target atoms to align with Ihe magnetic field. 
The scanner also includes coils tiiatqyply a transverse niagn RF pulses are 

emitted by the coils, causmg the target atoms to absorb energy. In response to the 
RF pulses, photons are emitted by the target atoms and detected as signals in . 

20 receiver coils. 

The signals detected in the recerver coils may be pxDcessed to construct 
image of the target The siggoals may be made proportional to flie spatial frequency 
content 0^-sp ace) of Ifae image throu|^ the proper qipfication of gradients to fibe 
magnetic field. Thek-spacemayconqpriseselsofsaniples,caUed l^^ 

25 corre^nding to a single phase 0icodingofthesanq;>ling process. SisweUknown 
tiiat the number and spacing of lines in k-space determines botii the field of view 
(^GV) and the spatial resohition of the reconstructed image. 

Data processing may be performed on the k-space sanq)les to produce a final 
image of die object in "^image space**, e.g. a spatial arrangement of pixels. The data 

30 processing is typically performed usmg a computer, which is any device conq)rising 
a processor and memory, wherein the processor executes instructions and acts upon 
data provided fix)m the memory. 
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R^id imagmg is desiiable in order to reduce the time lequiied to perform 
volume imagmg consisting of a large number of slices, to reduce the breath-hold 
time, or for dynamic imagmg applications such as fimctional imagmg of the heart or 
brauL R^id imaging also provides increased motion tolerance. A number of 
5 accelerated imaging methods have been developed. In several of these methods, 
undesirable "ghost! artifacts arise when the k-space samples are processed into the 
image domain. In one such method, known as echo-planar imaging (EPI), 
distortions in k-space roay lead to image domain ghosts. In another such method, 
known as SENSE, intentional k-space undersampling (sampling fewer lines than the 

10 number required to image a chosen field-of-view) accelerates the data acquisition 
but results in imag^ domain ghosts due to aliasing. For more details on SENSE see 
Pruessmann et al., SENSE: Sensitivity Encoding for Fast MRIy Magnetic Resonance 
in Medicme, 1999 Nov; 42(5): 952-962.) With the SENSE approadi, ^ort 
i?^ch arise fiom aliasing may be siqypiessed in the ima^ byway of a technique 

15 known as *^hased array combining^. Hbased array combining may be applied to 
siq;ypress ghost arti&cts arising fiom a variety of mechanisms^ not just aliasing. 
United States Patent Application No. 09/825,617, entitled G%ot^i47^^M 
Cancellation Using Phased Am^ Processing, and filed on .^iril 3, 2001, by 
KeUman et al. (hencefordi "KeUman V^^ teaches one such phased array combining 

20 qiproacli. 

Phased array ccmdnmng approadies for g^iost cancellation (mduding 
SENSE) involve combining multiple intermediate images, each comprising ghost 
arti&cis, to produce a fiboal image in wMch^iostarti£u4sax^ Often, the 

intermediate images are combined in a manner wfaidi is numerically ill-conditionfid, 

25 so that noise in the intermediate imag^ as well as errors in the combiidng wd^ts 
(brou^t on, for instance, by noise in the operation and characterization of the signal 
reception process) may aniplify noise in the final image. One technique to mitigate 
this problem is called regularization or matrix conditioning. Regularization involves 
a tradeoff between a level of ^ost arti&ct suppression and noise amplification. 

30 Current approaches s^ly a fixed amount of regularization to all pixels of an image. 
Such qiproaches do not take into account that particular regions (sets of one or more 
pixels) of flie final image may braefit fi:om substantially less g^st siq>pression than 



2 



•) 

wo 03/032816 



PCTAJS02/33571 



Others, and Aerefore can bmefit fiom greater noise reduction by trading offmore 
ghost siqipression than in other regions of the image. 

SUMMARY OF THE DISCLOSURE 

S In one aspect, a MRI imaging system inchides at least one processor and 

coils to acquire k-space samples of a target to image. The system includes a 
machine-readable media comprising instnictions which, when executed by the 
processor, result in determining different regularization matrices for different 
regions of an image of the target The regularization matrices are applied in the 

10 detenninationofunmixing matrices for the different pixel positioiis. The unmixing 
matrices are ^plied to produce the image without ghost arti&cts jBx>m intermediate 
MRI images produced from the plurality of k-sfpace sanq>les, the intennediate 
images each having ghost arti&cts. 

In another aspect, a MRI imagmg system includes coils to acquire reduced k- 

I S space samples for intermediate images. An apparatus comprising at least one 

processor executes instructions to combine the intermediate imaged produce a full 
field ofviewiniagp, by applying unmixing roatrices. Each immixing matrix is 
regularized according to a target level of alias suppression for a region of flie full 
field of view ima^. 

20 

BRIEF DESCimrnON OF THE DRAWINGS 

FIG. 1 isalug^levelblodcdiagramofanMRImettiodeinbodimeaa^ 
FIG. 2 is a more detailed blodk diagram of an MRI melbod embodiment 
HG. 3 is a block diagram showing an GDobodiinent of data processing to 
25 reduce ghost arti&cts in the final image. 

FIG. 4 is a block diagram showing an embodiment of an adaptive, spatially- 
variant regularization method. 

FIG. 5 is a block diagram of an ^aratus embodiment 

30 DETAILED DESCRIPTION 

In the following description and figures, like numbers refer to like elements. 
References to ''one embodimenf and ''an embodiment^ do not necessarily refer to 
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^ same embodiment, atthoii^ they may. Fi]i&emiare,1lie singular terms 
'^an," and "Hfae" include plural referents unless the context clearly indicates 
otherwise. 

Tbrou^out this description, embodiments may be described which apply 
5 matrix regularization on a pixel-by-pixel basis. Other embodiments may operate in 
substantially similar manners but may regularize according to regions of an image 
larger than a single pixel (for exdmplc, upon groups of neighboring pixels). In the 
figures, inoiages may be illustrated as pixel mqps. These pixel m^s are provided for 
illustration purposes only and may represent, in some instances, the complex values 
10 from which the pixel intensity magnitudes are generated. 

Overview 

HG. 1 is a hi^4evdblo€k diagram ofanMRImelhod embodiment 100. A 
set of k-£fpace samples is acquired firom Ne receiver coils, where Nc is an integor 

IS number. F<»r example, Nc =4 when four coils are onployed. Data processing is 
performed at 102 to produce a final image 104. In general, the data processing 
involves (1) producing a plurality of intetmedSate images finom the reduced k-^ace 
sanQ>les^ and (2) combining fiie intermediate images into a final image. The 
intennediate images may contain ghost acti&cts arising fixxm the k-q>ace acquisition, 

20 for exarnple arising fiomEPI distortions or undersanqpling. Hie intermediate . 
images maybe conibined to produce the final image ^vfdnle siqq;nnessnig ^bost 
arti&cts. 

For simplicity, the following discussion appUes to an application conq[msing 
a smgle^ostarti&ct in each intermediate image. The discussion is gmerally 
25 qyphcable to situations involving a greater number of gjbiost arti&cts as well. 

Acquiring the Intermediate Images 

FIG. 2 is a more detailed block diagram of a MRI method embodiment 200. 
MR coils (in this exanq)le, coill» coil2, coil3, and coil4) may be positioned in 
30 differing aspects relative to an object to sample. K-space samples may be acquired 
for each coil. At 201 , 202, 203 , and 204, an EFI operation is ^lied to each set of 
acquired k-space samples, resulting in intermediate images 205, 206, 207, and 208, 



4 



wo 03/032816 



PCT/US02/33571 



respectively* Id this example, each of the intetmediate images 205-208 is compiised 
of ttie desiied image plus a single ghost artifact The ghost arti&ct in this example is 
a shifted version of the desired final image, with a vertical shift equal to FOV/2 
pixels. This situation may arise in SENSE applications where image acquisition is 

5 accelerated by a factor R=2, where each set of k-space samples comprises FOV/2 
lines, as well as in other applications such as EPL The intermediate images 205-208 
in this example con^>rise pixels which are a siqjeiposition of a desired pixel value 
and a ghost pixel value. The number of superpositions at a pixel location is equal to 
R, the acceleration factor. Thus, for R=2, a pixel position (x,y) in eadi mtermediate 

1 0 image noiay comprise pixel values which are a superposition (sum) of (1) the pixel 
value at position (x,y) in the final image, and (2) the pixel value &om a position 
(x,y+-FOV/2). (More specifically, position (x, Cy+FOV/2) mod FOV) where FOV is 
the pixel hd^ of flie image. Henceforth, the position (y+FOV/2) mod FOV) will 
beiefeuedtosin^lyasyaiiasO This ghosting efifect is also lefenred to as "^Id over" 

15 and "tnixing^. 

Data processing 1 02 may be s^lied to combine the intermediate images 
205-208 to produce the final image 104. When cornbining, proper weighting of fiie 
pixel values of the intermediate images 205-208 may produce pixels in the final 
image 1 04 with suppressed ghost artifacts. 

20 

Coiiibining the Intermediate Images 

FIG. 3 is a block diagram showing one embodiment 300 of data processing 
to siq[>pressg^startifiicts in ttie find image. Pixel values pi(x»y)finom 
corresponding positions (x,y) of the four intermediate images 205-208 maybe 

25 processed to produce a final image 104 with suppressed ghost artifacts^ where i has 
therange 1 to4(inthisexanq;>lewith4coils). Thn)ugh application of an unmixing 
matrix at 306, fiie pixels pi(x,y) of the intermediate images 205-208 may be 
wei^ted and added together to produce a pixel q(x,y) of the final image 104. Noise 
may be introduced into the final image due to both (1) noise in the received k-space 

30 samples, and (2) inaccuracies in the weights sqipUed to the pixels of the intermediate 
images. The iuaccuiades in the weights, in turn^inay result fiomiriaccurad 
estimated MR coil sensitivities, among other things. The level of noise introduced 
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into the £uialixnage is dependent upon tt^ The noise may be reduced at the 

expense of reduced artifact suppression. 

For Nc=^, one pixel fix>m each position (x,y) of the intennediate images 205- 
208 are arranged into a 4x1 (Ncxl) vector P. .A|)plication of the unmixing matrix 
5 may transform P into matrix Q comprising pixel vahies of the final image 104. In 
other words, the intermediate images 205-208, when taken together, comprise 
enough information to separately resolve all pixel values in the final image 104. 

The transformation of P to Q, e.g. the separation of superimposed pixels of 
tibie int^mediate images 205-208 into non-siqperinq)osed pixels of the final imagq 
10 104» may be accomplished by q>plication of the unmixing matrix U. 

Q=U*P 

hi the SENSE accelerated imaging q^plication, ^eare intermediate images are 

15 comprised of R siq>erinq>osed aliased pixels, the unmixing matrix U is typically 
formulated as a 2xNc (RxNc) with pixels locations (x,y) evaluated over FOV/2, 
which produces a2xl (Rxl) vector Q, with elranents correqxmding to (x,y) and 
(x^yaiias). The unnaixingmatrixmayequivalenUybe formulated withU defined a^ 
IxNc vector and pixels locations (x»y) taken Scam fbe fidl FOV. 

20 A result of the matrix multiplication of U and P is that eadh demmt of Q 

represents a wdgbted sum of Nc pixels, one fiom each corresponding position (x,y) 
ofthe intermediate images 205-208. The elements ofU are detennined such that 
fids weighted sum suppresses ^sting effects which may be present in the pixels of 
the intennediate images 205-208. hi other words, fiie weighted sum produced by U 

25 for each pixel position (x,y) of the intennediate images 205-208 sq>arates the pixel 
values at that position. The degree ofsuch separation may vary, &r reasons and in 
manners described below. 

The elements of U may be detennined &om transformations involving an 
NcXRmatrixS. Matrix S is referred to as the "sensitivity matrix". EachrowiofS 

30 corresponds to the "sensitivity" ofcoil lata plurality of pixel positions. ForR==2, 
each column j of S comprises (1) a coil sensitivity value for the pixel at position 
(x,y) of the final image 104, and (2) a coil sensitivity value fi>r the pixel at position 
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(x,ya}ias) of the final image 104. The pixels at positions (x,y) and (x, yaiias) of the 
fiaal image are the pixels which superimpose at position (x,y) in the intennediate 
images 205-208. 

Assuming that the elements of S perfectly represent of the actual sensitivities 
5 of the coils, then one manner of determining the unmixing matrix U from S involves 
applying the transfonnation 

10 where rq)resents the conjugate transpose (Hermitian operation) of matrix S» the 
negative exponent (- 1) indicates a matrix inversion oporation^ and iir represents the 
well-known noise covariance matrix for tibe receive coils. Here, Uo represents an 
intennediate determination of U which may comprise an undeshable gain aspect. 
Manneis of compensating for diis gain aspect to produce the final matrix U are more 

15 fiiUy described below. 

Typically, tiie elements of S are not perfect rqsresCTtations of the actual 
sensitivities of the coils. It may not be possible to experimentally detmnine tiie coil 
sensitivities with complete accuracy for a particular imaging operation. 
FurtfaCTnore^ the sensitivity values may vary according to interactive effects 

20 between coils, motion ofthe target object, and other variables. Thus, tilie matrix S 
may comprise errors which deviate fitmx the actual soisitivity values ofthe coils. 
Furthermore, the k-space samples may comprise noise which may add to die pixel 
values ofthe intermediate images 205-208. Application ofthe unmixing matrix U to 
the pixels ofthe intermediate images 205-208 may amplify this noise, and enx>rs 

25 resulting fixnn the errors in S, to undesirable levels in the final image 104. 

It may be possible to reduce some of die effects of noise, and enors in S, by 
applying regularization to the detenniiiation of U, in manners which are further 
described below. Furthermore, an ads^tive approach to determination ofthe 
sensitivity values may lead to reduced errors in S and thus improve ghost artifiict 

30 suppression in the final image 104. 
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Determination of Coil Sensitivity Values 

Id one embodimeoot, coil sensitivity values may be detennined by acquiring 
full FOV reference images for each coil. These reference images may be free of 
substantial ghost artifects. These reference images may be acquired either before or 
5 after the desired imaging operation (for example, accelerated imaging) by 
performing a 'reference scan'. 

In another embodiment, the reference images noiay be acquired, at least in 
part, during the course of the imag^ig operation. One such qiproach is described in 
U.S. Patent Application No. 09/735,263, entitled^cce/era/^/ Magnetic Resonance 

10 Imaging, and filed on December 1 1, 2000, by Kellman et al (henceforth ''KeUman 
2*0. Kelhnan 2 teaches a niamier in ivfaidifuU FOV k-space sample 
acquired during the course of an accel^ated imaging operation and processed into 
reference images. The reference images may ften bellied to adq)tively 
detennine the sensitivity values of tibe coils at different times during the course of 

IS the accelerated imaging operation, hi other words, botii accelerated reduced k-space 
sanq^hng and slower, fiill FOV k-space sampling take place over the course of the 
imaging operation. The iniages derived fiom the fiiU FOV samples are q[>pUed to 
ad^tively detennine sensitivity values for the recdver coils. 

In this maimer, imaga tar^ motion and other environmental dhanges that 

20 may effect coil sensitivities maybe taken into account during flie course of 

acceletated inu^ED^ The reference ima^aie Ml FOV, and thus take longer to 
acquire, flian do Ihe accelerated reduced k-space images, hi other words, tfie 
reference images have a lower temporal tesohftion than reduced k-space images, and 
among other distortions may contain blurring effects due to is^d image target 

25 motion (for example, the ventricular motion ofa beating heart). However, coil 

sensitivities adaptively estimated fi^om such reference images may nonettieless prove 
more accurate over the course of the imagmg operation than sensitivities estimated 
once, prior to the imaging operation. 

hi addition to ad^tively estimating the coil sensitivities, it may be possible 

30 to further reduce noise effects in the final image 1 04 by cornpromising some ghost 
artifiict suppressiocL Noise effects in the filial image 104 may be further reduced, at 
the expense of ghost arti&ct siqipression, by "regularizing^ or ^^better conditioning^ 
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the matrix inverse operation (S^'^S)'^ in the determination of matrix U. One 
manner of regularization involves adjustments to the elements of the term (S^'^S), 
such that inversion of the temi results in less noise amplification. Using this 
approach, the regularization may be performed by the addition of a matrix A» taking 
S the form 

SV^S+A 

Thus, the regularized determination of matrix U becomes 

10 

Uo = (SV*S+A)"^S"\|;"^ 

Again, Uo represents a determination of U whidbi may comprise an undesirable gain 
asq[>ectandwMchmaybeconipensatedforinmamier5tobedescxi^^ Priorait 

lis regularization techniques have q>plied a diagonal RxR matrix A in which the 
diagonal elements of A have constant values near but slightly greater tlian the 
smallest eigenvalues of the term S^"^S. See for example SENSE Image Quality 
Improvement Using Matrix Regularization, KLF. King et al-. Proceedings of the 
International Society of Magnetic Resonance in Medicine 9, 1771 (2001). In the 

20 prior art, the elements of matrix A are constant for all pixel positions, and are not 
adsQitively detemxmed during the course of an iinagjm 

Adaptive Deterininatioii of Rq^nlarizatioii Values 
FIG. 4 is a block diagram showing an embodiment 400 of an adiq[>tive, 
25 spatially-variant r^;i]larization method laoneeoibodiment, AisadiagonalRxR 
matrix with element values detennined accordmg to a targ^ level of g^ost arti&ct 
suppression to apply to one or more pixels in the final image. The target level of 
ghost arti&ct siQ>pression to ^ly, and the coil sensitivities, may be ad^tively 
determined according to reference images acquired during the course of an imaging 
30 operation, as weU as static reference images acquired before and/or after 
operation. 
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lot one embodiment, a set of reference images ij (one for each coil fiom 
which k-saiiq>les are acquired) is applied to determine fhe sensitivity matrix S at 
404. The reference images ri may be combined to produce a *combined-m^nitude' 
reference image rem which is ^plied to detennine, at 402, a target level of £^05t 

5 artifact suppression to ^ly to each pixel of the final image. In one embodiment, 
the value of a pixel rcm(x,y) comprises the combined magnitude of tiie corresponding 
pixel values in ri, e.g. the square root of the sum of the squares of the complex pixel 
values ri(x,y). In another embodiment, for which generalized phased array ghost 
cancellation produces R separated ghosts (as described in Kellmanl), the combined 

10 magnitude for each individual separated ^ost image may be input to the 

detenninationoftarget ghost suppressioii, at 402. In this case, the target arti&ct 
sq>pression may be calculated as the ratio of desired pixel and ghost pixel using 
pixels in separated images. The sensitivity matrix S maybe calculated from either a 
single image (as shown in FIG. 4) or fiom multiple separated ghosts as described in 

IS Kellmanl . In another embodiment, the combined magnitude lefercnce image 

loi^^y) used for calculating the regularization matrices may be calculated fiom an 
initial application of unmixing using unmixing matrices (Uo) calculated using a 
snoiaUer fixed regularization or without any regularization. 

Consider a pixel at position (x,y)ofthe final image. ForR=2,the 

20 corresponding pixel in the intennediate bnagBs will have a value viAidi is file 
superposition (sum) of the pixel values in fiie final image at positions (x,y) and 
(x,yaiias)- Thus, in one embodnnent, file target level of ^ostartifitt*siq)pressiOT 
^ly to file pixel at position (x,y) of tiiie final image is proportional to the ratio 
rcn,(x, yaKa8)/rcm(x,y). For exan?)le, if fiie pixel at position (x,y) has an intensity vahie 

25 of one (1), and the intensity value of the pixel at position (x, yaKas) is five (5), the 
ratio is5:l. A target level of ghost artifact siq?pression to ^ly to the pixel at 
position (x,y) may be around five times larger than a target level to ^ly to a pixel 
at a position whore the ratio is closer to one. 

The taiget levels of ghost arti&ct suppression may be chosen according to 

30 various criteria, to one anbodtment,flie target levd of gjiost artifact suppression to 
apply to a pixel at a position (x,y) is chosm as a percentage of the intensity value of 
tiie pixel in the refermce image. In another embodiment, a level of noise present in 
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a pixel value at position (x,y) of the ref^pence image is determined. A target level of 
ghost artifact suppression is chosm to reduce tiie g^ost artifact of the pixel to a level 
on order with the level of noise. 

Ctoil sensitivities are estinaated at 404. The coil sensitivities may be 

5 estimated from the reference image in various mamiers, for example in the manners 
detailed in KeUman 2. The estimated coil sensitivities S are provided to determine 
the unmixing matrix U at 408. 

At 406, the regularization matrix is deteraiined according to tiie target levels 
of ghost arti£ict suppression. In one embodiment, initial values are chosen for the 

10 elements oflbe regularization matrix A. At 410, actual levels of ghost arti&ct 
siqjpression for each pixel of the final image axe determined fiom the matrix U 
which was detennined according to (l)flus initial A, and (2) the estm The 
noise covariance maybe incorporated for optimized SNR as in previous described 
matrix formulation for U. The noise covariance may be a separate noise-only 

IS reference scan or estimated daring imaging fiom noise-only pixels. Iterative 

adjustments to the elements of A are made until the target level of ^ost artifact , 
suppression is achieved For each iteration, fiie actual level of ^bost artifact 
siqppression is checked with the target level of ghost arti&ct suppression for the 
pixel, and the element values of A are adjusted accordingly, if necessary, to bring 

20 the actual level of ghost arti&ct siq>pression closer to the target level. 

An actual level of ghost arti&ct suppression for a pixel at position (x,y) of 
the final image may be determined &om the matrix product, 

p=U#S 

25 

^a:e tihe diagonal of p may be nonnalized to all ones (Is). For R=2,fte off- 
diagonal elements of p are thGti each proportional to the actual ghost arti&ct 
suppression ^Ued to fiie pixels at positions (x,y) and (x, yaiias). For R=^ 
adjustments to the diagonal elements of A may indq>ende(ntly affect the elements of 
30 p. In other words^ adjusting a particular diagonal element of A niayindq)^ 

affect the actual ghost arti&ct suppression for a particular pixel of the final image. 
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For R=^ adjusting element Ki^ affects the actual g^ost artifiict suppression 
represented by element pi^, independent of p2,i. Likewise, adjusting element Ai,i 
affects the actual gliost artifact suppression represented by element p2,i, independent 
ofpi;2. 

In one embodiment, this process may be repeated for each pixel of the final 
image until &e actual ghost arti&ct suppression is in accordance wiOi the target 
level. Qfcourse, other numnersofdetenzuning the values of A whi(Aa^ 
desired ghost aitiiGu^t suppression may also be employed, such as closed-end 
solutions, which do not involve an iterative process. 

At various stages of tihe imagjxig operation, a new reference image may be 
provided to iq>date the determination of target ghost aiti&ct suppression levels and 
coil sensitivities. See Kelfanan 2 for a description of one manner in which this may 
be done. To compensate for motion distortions in the low temporal resolution 
refisrence images, an ^order filter' may be applied in one embodimenL 
Corresponding elements of tfie determined regularization matrix A may be coiqpared 
within a neighborhood of pixels. The elements of the A matrix for the pixel at the 
center of the neigbbodiood may then be set to the minimum element values of the A 
matrices for all pixels in the nei^bLboifaood. For exanq[>le, consider the following 3x3 
pixel neagbbodiood 

PlP2p3 
P4PSP6 

The matrix A for ps may be set to coiiqirise the niinmmm element values of A for 
the set of pixels pi-p9. This process may be repeated for each pixel of the final 
image to reduce errors related to the low tenq[K)ral resolution of the reference image. 

Adjusting the Gain 

RecaU tiiat in one em]x>diment, an umnixing rnatrix Uo coin^^ 
zsfGcX is determined by 
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Uo=(sVs+A)-^sV 

In one embodiment, a gain matrix G may be included into the detennination of U> to 
S conopensate for liie gain aspect of Uo, as follows, 

U = G*Uo = G(S VS + A)'^S V* 

For a diagonal matrix A, the terai G(S^"^S + A)"* £q[»proaches G(A)'* as the values 
10 ofthe diagonal elements of A are increased. Inoneembodime3]t»thetennG(S^~^S 
+ Ay ^ be reduced to G(Ay^ when the diagonal values of A are substantially 
larger than the maximum eigenvalues of tiie temi S^'^S. In this case, 

U = G(A)"^S V* (approximately) 

15 . 

It is known that, in tiie absence of ghost arti&cts, optimal signal-to-noise ratio 
(SNR) may be achieved in the final image where U is proportional to ^ 
However, in this case no ghost artifact suppression takes place. See for exanq)Ie The 
NMR Phased Array, Roemer et aL, Magnetic Resonance in Medicine 1990; 16:192- 
20 225 (henceforth Roemer). Li one embodiment, tiie target ghost artifact 
sqypression for apixel is low» tiie tean G{A)'^ maybe set to ^^proximately the 
identity matrix I, in which the vahie of all diagonal elements is approximate!/ one. 
Thus, 

25 u=G(Ay^sV=isV^V^ 

In otiier words, for pixels where ghost arti&ct siqipression is not substantially 
needed, tiie elonmts of tiie matrices G and A may be determined such that the SNR 
for the pixel in the final image is close to the optimum levels as determined, for 
30 exanqple, in Roemer. 
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For pixels where the ghost ariifict suppression to apply is substantial, G may 
be detennined such that Hxe diagonal elements of tib^e matrix product G«p are close to 
one. In other words, 

J- 0 

Pll 

0 -L 

P2IJ 

5 where pn and P22 are diagonal elements of the matrix po= Uo»S. As a consequence 
of compensating for the gain Bspcct of Uo, the pixels in the final image may have a 
more desurable intensity. 

FIG. S is a block diagram of an qypaiatus mbodiment 500. The sqpparatus 
500 comprises a processing miit 502 (e.g., a processor, microprocessor, micro- 
10 controUer, etc.) and machine-readable media 504. Depending on the configuration 
and qjplication (mobile, desktop, server, dtc.X the memory 504 maybe volatile 
(such as RAM), non-volatile (such as ROM, flash memory, etc.) or some 
ccnnbinatton of the two. By way of sample, and not limitation, machine readable 
media 504 may conqirise volatile and/or nonvolatile media, removable and/or non- 
15 removable media, including: RAM, ROM, EEPROM, flash memory or other 
memory technology, CD-ROM^ digital versatQe disks (DVD) or other cp&cal 
storage, magnetic cassettes, magnetic tape^ magnetic disk stora^ or otfao: magnetic 
storage devices, or any oHker medium wbich can be used to store tte desired 
infinmation to be accessed by the apparatus 500. The xnadbdiie readable media 504 
20 may be ixapkmessted in any mediod or technology for storage of information such as 
conq)uter readable instructions, data structures, program modules or other data. 

The media 504 may comprise instructions and/or data which, when executed 
by the processor 502, may result in the apparatus 500 carrying out acts in accordance 
with embodiments of the methods described herein. 
25 The ^q>aratus 500 may comprise additional storage (removable 506 and/or 

non-removable 507) such as magnetic or optical disks or t^. The qrpaiatus 500 
may further com^nise input devices 510 such as a keyboard, pointing device, 
microphone, etc., and/or output devices 512 such as display, speaker, and printer. 
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The apparatus 500 may also typically include netwoik coimectiotis 520 (such as a 
uetwoik ad^ter) for coiq)liiig to other devices, cooaputers, networks, servers, etc. 
Using eitiier wired or wireless signaling media. 

The coiuponents of the device may be embodied in a distributed conq>uting 
5 system. For example, a temiinal device may incorporate input and output devices to 
present only the user interfece, whereas processing component of tiie system are 
resident elsewhere. Likewise^ processing functionality may be distributed across a 
plurality of processors. 

The ^aratus may generate and receive machine readable instructions, data 

10 structures, program modules or other data in a modulated data signal such as a 
canier wave or other transport mechamsoL These instructions and/or data may, 
when executed by tfie processor 502, result in acts in accordance with procedures of 
file present invention. The term **modulaled data signal" means a ^gnal that has one 
or more of its chaoractedstics set or changed in such a manner as to encode 

15 information in the signal Byway of exanoplei, and not liimtati(m,commuiiication 
media includes wired media such as a wired netwoik or direct-wired coimection, and 
wireless media such as acoustic, RF, infrared and oflier wireless media. 
Conununications media, including combinations of any of the above, should be 
understood as within the scope of machine readable media. 

20 Having described and iUnstrated the present invention with refermce to one 

or more illustrated embodiments, it will be recogjtuzed that the illustrated 
embodiments can be modified in arrangement and detail without departing &am the 
principles and scope of fiie preseat invention. It should be understood that ttie 
programs, processes, or methods described herein are not related or limited to any 

25 particular type ofdata processing ^yparatus, unless indicated othermse. Various 

types of general purpose or specialized data processing q>paratus, including desktop 
computers and workstations, may be used with or perform operations in accordance 
with the teachings described herein. Elements of the illustrated embodiments may 
be implemented in software, hardware, firmware, or combinations thereof 

30 In view of the many possible embodiments to which the principles of the 

present inv^tion may be apphsd^ it should be recognized that the detailed 
embodiments are illustrative only and should not be taken as limiting in scope. 
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Ridher, the present inventioii eacompasses all such embodiments as may come 
within the scope and spirit of the following claims and equivalents thereto. 
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WeQaim: 

1. A compiiter-m^lemeiited inefhod, conqmsmg: 

detemiining a plurality of different regulaiization matrices for a plurality of 
difTerent regions of an image; 

flying the regularization matrices in the detennination of a plurality of 
unmixing matrices for the differait regions; and 

flying the unmixing matrices to a plurality of undersampled MRI imag^ 
with g^ost arti&cts to generate the image without ghost arti&cts. 

2. The method of claim 1 wherein d^enniniiig a plurality of different 
regularization matrices furtfa^ comprises: 

ad^tively determining the different regularization matrices according to a 
plurality of reference images acquired over the course of an imagmg operation. 

3. The mediod of claim 2 whearein detennining a plurality of different 
regularization matrices further comprises: 

detecmining the regularization matrices according to a plurality of different 
target g^ost siqypression levels for the different regions of the imagp. 

4. The mediod of claim 1 wherefai the detenniiiatiQn of a plurality of 
unmixing matrices for tiie plurality of different regions further conprises: 

adaptively detenninrng coil sensitivhy values according to a plurality of 
reference images acquired over the course of an imaging operation; and 

qyplying the adaptively determined coil sensitivity values to the 
determination of the urmiixing matrices. 

5. The method of claim 3 wherein determining a plurality of different 
regularization matrices further comprises: 

ad^ting the pluraUty of regularization matrices according to differences 
between actual ghost siq>pression levels for the plurality of pixel positions and the 
plurality of different target gjiost sq>pression levels for the plurality of regions. 
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6. The metliod of claim 1 wherein determiniiig aphitality of different 
regolarizatioii matrices far a plurality of different regions of an image &rther 
comprises: 

detennining the regulaiization matrix for a first region of the different 
regions fiom the minimum element values of the regulaiization matrices of 
surrounding regions. 

7. The method of claim 1 in which each region comprises a single pixel 
of the image. 

8. The method of claim 6 in whidi eadi region comprises a single pixel 
of the image. 

9. The mdhod of claim 3 in wfaidi the plurality of taig^ g^ost 
suppression levels are chosen as a percentage of an intensity value of the different 
regions. 

10. Themelhodof claims in which a level of noise present in the 
different regions is d^ermined and the plurality of target gjbost suppression levels 
are chosen to reduce a g^ost artifict of the different regions to a level on order with 
die level of noise. 

11. A method conqnising: 

acquiring a plurality of undersairqpled MRI images; and 
combining the pluraUty of undersanopled MRI iiimges to produce a full field 
of view image by applying a plurality of mmaixing matrices, each unmixing matrix 
regularized according to a target level of alias siq>pression for a region of the fiill 
field of view image. 

12. The me&od of claim 1 Ifurther conipiising: 

adjusting tfie target level of alias si^>pression accordmg to reference unages 
acquired over the course of an imaging operation. 
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13. An MRI imaging system comprising: 
at least one processor, 

a plurality of coils to acquire a plurality of k-space samples of a target to 
S image; and 

a machine-readable media conq>rising instructions whidi, when executed by 
the processor, result in 

determining a pltirality of different regularization matrices for a plurality of 
different regions of an image of the target; 
10 applying the regularization matrices in the determination of a plurality of 

unmixing matrices for the regions; and 

applying tiie umnixiag matrices to produce the image without g^ost arti&cts 
fiom a plurality of MRI images produced fiom flie plurality of k-space samples and 
each coniprisiiig gjiost arti&cts. 

15 

14. The system of claim 13 wherein the iostractions, when executed by 
the processor, further result in: 

ad^tivdy determining the diff^arent regularization matrices according to a 
plurality of reference images acquired over the course of an imaging operation of tibie 
20 target 

15. The system of claim 14 wh&xm the instructions, vfbm executed by 
the processor, iiir&er result in: 

determining the regularization matrices according to a plurality of different 
25 target ^ost siqypression levels for the regions of the image. 

16. The system of claim 13 wherein the instructions, when executed by 
the processor, furtibter result in: 

adaptively deterinining sensitivity values for the coils according to a plurality 
30 of reference images acquired ovct the course of an imagmg operation of Ifae target; 
and 



19 



wo 03/032816 



PCTAJS02/33571 



q[yplymg the adq>tively detennined coil sensitivity values to the 
detennination of the unmixiag matrices. 

17. The system of claim 15 wherein the instructions, when executed by 
the processor, further result in: 

adapting the regularization matrices according to dijSerences between actual 
^ost suppression levels for the regions and the different target ghost suppression 
levels for the regions. 

18. The system of claim 13 wherein the instmctions, when executed by 
the processor, further result in: 

detemuning the regularization n:iatrix for a first region of the difTermt 
^regions fiom the minimum element values of the regularization matrices of 
sunounding regions. 

19. The system of claim 13 wherein eadi different region comprises a 
single pixel position of the image. 

20. The system of claim 18 wherein each different region comprises a 
singlle pixel position of the image. 

21. An MRI imaging system, coniprising: 

a plurality of coils to acquire reduced k-space sanqples for apturality of 
intermediate images; and 

an apparatus conq)rising at least one processor to «ecute instructions to 
combine &e intermediate images to produce a full field of view image by q^iplyiog a 
plurality of unmixing matrices, each unmixing matrix regularized according to a 
target level of alias suppression for a region of the fiill field of view image. 
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22. The system of claim 2 1 wherein fhe coils fiiither operate to acquire at 
least one reference image during the course of an imaging operation, and wherein 
the instructions, vifbesi executed, fiirtherresutt in adjusting the target level of alias 
5 siqppression according to the at least one reference images. 
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